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ABSTRACT 

An analysis of the environments around a sample of 28 SCR radio galaxies with red- 
shifts 0.6 < z < 1.8 is presented, based primarily upon K-band images down to 
K ^ 2{) taken using the UK Infrared Telescope (UKIRT). A net overdensity of K- 
band galaxies is found in the fields of the radio galaxies, with the mean excess counts 
being comparable to that expected for clusters of Abell Class richness. A sharp 
peak is found in the angular cross-correlation amplitude centred on the radio galax- 
ies, which, for reasonable assumptions about the luminosity function of the galaxies, 
corresponds to a spatial cross-correlation amplitude between those determined for low 
redshift Abell Class and Abell Class 1 clusters. 

These data are complimented by J-band images also from UKIRT, and by optical im- 
ages from the Hubble Space Telescope. The fields of the lower redshift {z < 0.9) radio 
galaxies in the sample generally show well-defined near-infrared colour-magnitude 
relations with little scatter, indicating a significant number of galaxies at the redshift 
of the radio galaxy; the relations involving colours shortward of the 4000A break show 
considerably greater scatter, suggesting that many of the cluster galaxies have low 
levels of recent or on-going star formation. At higher redshifts the colour-magnitude 
sequences are less prominent due to the increased field galaxy contribution at faint 
magnitudes, but there is a statistical excess of galaxies with the very red infrared 
colours {J — K > 1.75) expected of old cluster galaxies at these redshifts. 
Although these results are appropriate for the mean of all of the radio galaxy fields, 
there exist large field-to-field variations in the richness of the environments. Many, 
but certainly not all, powerful z 1 radio galaxies lie in (proto-)cluster environments. 

Key words: Galaxies: clustering — Galaxies: active 



1 INTRODUCTION 

Clusters of galaxies are the largest, most massive, collapsed 
structures in the Universe, and as such are of fundamental 
importance for many cosmological studies. They provide a 
unique probe of large-scale structure in the early Universe 
and, as the systems which separated from the Hubble flow at 
the earliest epochs, they contain the oldest galaxies known. 
These can strongly constrain the first epoch of the formation 
of ellipticals, and hence set a lower limit to the age of the 
Universe; the effectiveness of using such old galaxies at high 
redshifts towa rds th is goal has been well demonstrated by 
Dunlop et al. (1996). Further, because clusters contain large 
numbers of galaxies at the same distance, they are important 
testbeds for models of galaxy evolution. 
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The cores of optically-selected clusters are dominated 
by a population of luminous early-type red galaxies which 
occupy a narrow locus in colour- magn itude relations. Stan- 



1998) showed that, out to 
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ford, Eisenhardt and Dickinson (1 
redshifts 2 ~ 1, the evolution of the colours of early-type 
cluster galaxies on these relations is completely consistent 
with passive evolution of an old stellar population formed 
at high redshift. The small intrinsic scatter of the galaxy 
colours, and partic ularly the fact that i t remains small at 
redshifts z > 0.5 (Stanford et al. 199S), implies that the 
star formation of the ellipticals comprising a cluster must 
have been well synchronized, and sets tight limits on the 
amount of recent star formation that might have occurred 
(e.g. Bower, Kodama and Terlevich 1998). 

Cluster ellipticals also show a tight relationship between 
their effective radius, effective surface brightness, and cen- 
tral velocity dispersion (the 'fundamental plane'; c.f Dressier 
et al. 1987, Djorgovski and Davies 1987). The location 
of these ellipticals within the fundamental plane has been 
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shown to evolve with redshift out to z = 0.83, in a manner 



dista nt radio sources, selected by infrared colour (Dickinson 



which implies that the mass-to-light ratio of the galaxies 



1997 



evolves as Alog(Af/I/r) ~ — 0.4Az, roughly in accordance 
with passive evolution predictions (e.g. van Dokkum et al 
1998 and references therein). 

These results are in qualitative agreement with 'mono- 
lithic collapse' models of galaxy formation, in which an el- 
liptical galaxy forms the majority of its stars in a single 
short burst of star formation at an early cosmic epoch. On 
the other hand, the hierarchical galaxy formation models 
favoured by cold dark matter cosmologies (e.g. White and 
Frenk 1991) predict a later formation of galaxies with star 
formation on-going to lower redshifts. These are supported 
by the appearance of a population of bluer galaxies in many 
clusters at redshifts z > 0.3 (the Butcher-Oemler effect; 
Butcher and Oemler 1978) and by the high fraction of merg- 
ing gal axies seen in the z — 0.83 cluster MS1054— 03 (van 
Dokkuii et al. 1999). Since finding clusters at high redshift 
is strongly biased towards the very richest environments, 
in which early-type galaxies will have formed at the very 
highest redshifts, hierarchical models can still explain the 
apparent passive evolution and small scatter of the colour- 
magnitude relation out to redshifts 2; ~ 1. 

It is clearly important to extend cluster studies out to 
still higher redshifts, but the difficulty lies in the detection 
of clusters. At optical wavelengths, the contrast of a cluster 
above the background counts is minimal at these redshifts: 
the deep wide-area ESO Imaging Survey (EIS) has found 
12 'good' cluste r candidates above re dshift 0.8, but none 
above z = 1.2 (Scodcggio et al. 1999). Addition of near- 
infrared wavebands helps (e.g. Stanford et al. 1997), but is 
still a relatively inefficient method. Selection using X-ray 
techniques is more reliable, but X-ray surveys are currently 
sensitivity limited: the ROSAT Deep Cluster Survey found 
about 30 clusters with 2 > 0.5, but none above redshift 



0.9 (Rosati et al. 1998). Chandra and XMM will make a 



big improvement here. An alternative approach is to use 
powerful radio galaxies as probes of distant clusters: these 
can be easily observed out to the highest redshifts, and there 
is growing evidence that at high redshifts they lie in rich 
environments. 

At low redshifts powerful double radio sources (FRII's; 
Fanaroff & Riley 1974) are associated with giant ellipti- 
cal galaxies that are typically the dominant members of 
galaxy groups; the only nearby radio source of compara- 
ble radio luminosity to the powerful high redshift radio 
gala xies is Cygnus A , and this source lies in a rich clus- 
ter (Owen et al. 1997). At a redshift z ^ 0.5, analysis of the 
galaxy cross-correlation function around FRII radio galax- 
ies (Yates, Mi ller and Peacock 1 989) and an Abell clustering 
classification (Hill & Lilly 1991) have shown that about 40% 
of radio sources are located in clusters of Abell richness class 
or greater. At 2 ~ 1, the circumstantial evidence that at 
least some powerful radio sources are located at the centres 
of clusters is overwhelming and includes the following: 



• Direct detections of companion galaxies with narrow- 
band imaging (McCarthy, Spinrad and va n Breugel 1995) 
and spectroscopic studies (Dickinson 1997). 

• The observation that powerful radio galaxies are as lu- 
minous as brightest cluster galaxies at z ~ 0.8, and have ra- 
dial light profiles which are well-matched by de Vaucouleurs 
law with large (10 to 15kpc) characteristic radii (Best, Lon- 
gair and Rottgering 1998b). 

• The radio sources display large F araday depolarisation 
and rotati on me asures (e.g. Best et al. (1998a); see also Car- 
illi et al. (1997) for redshift z > 2 radio sources), requiring 
a dense, ionised surrounding medium. 

• Theoretical arguments that to produce such luminous 
radio sources requires not only a high AGN power, due to 
a very massive cent ral black hole being fueled at close to 
the Eddington limit (Rawlings & Saunders 1991), but also a 
dense environment to confine the radio lobes and convert the 
jet kinetic energy efficiently into radiation (e.g. see Barthel 
and Arnoud 1996). 

At still higher redshifts, radio sources have bee n de- 
tected out to redshift z = 5.2 (van Breugel et al. 1999), and 



some well-studied sources are known to lie in cluster envi- 
ronments. For example, towards the radio galaxy 1138—215 
(z = 2.2) extended X-ray emission has been detected (Car- 
illi et al. 1998) and n arrow-band ima, ging reveals over 30 
nearby Ly-a emitters (Kurk et al. 2000). Radio sources may 
therefore offer a unique opportunity to study dense environ- 
ments back to the earliest cosmic epochs. As yet, however, 
there have been no systematic studies of radio galaxy envi- 
ronments much beyond z ~ 0.5, and so it is important to 
investigate in detail the nature of the environments of the 
general population: do all powerful distant radio galaxies 
lie in cluster environments, or do only a minority but these 
have 'grabbed the headlines'?; what are the properties of 
any clustering environments (richness, radius, shape, etc.) 
surrounding these objects?; what is the nature of the con- 
stituent galaxies (morphological composition, segregation, 
etc.) of any detected cluster? 

In this paper, deep near-infrared observations of the 
fields of a sample of 28 powerful radio galaxies with 0.6 < 
z < 1.8 are analysed, in conjunction with Hubble Space 
Telescope (HST) images of the same fields, to investigate 
the ubiquity and richness of the environments of the ra- 
dio galaxies. The observations, data reduction, and source 
extraction are described in Section ^ In Section the in- 
tegrated galaxy counts are considered, in Section the an- 
gular and spatial cross-correlation amplitudes are derived, 
and in Section |5| an investigation of the colour-magnitude 
and colour-colour relations is carried out. The implications 
of the results are discussed in Section ^. Throughout the 
paper values for the cosmological parameters of f2 = 1 and 
Ho = 50kms~^Mpc~^ are assumed. 



• Detections of luminous X-ray emission from the fields 
of the radio galaxies (e.g. Crawford and Fabian 1996), some- 
times observed to be extended, indicating the presence of a 
relatively dense intracluster medium. 

• Large over-densities of galaxies in the fields of some 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 The dataset 

The data used for this research wa s pres ented and described 
by Best, Longair and Rottgering ( 1997 ). In short, the sam- 
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pie consists of 28 radio galaxies with redshifts 0.6 < z < 1.8 
drawn fr om t he revised SCR sample of Laing, Riley and 
Longair (198J). The fields of these radio galaxies were ob- 
served at optical wavelengths using the Wide-Field Plane- 
tary Camera II (WFPC2) on the HST generally for one orbit 
in each of two different wavebands. They were also observed 
at near-infrared wavelengths using IRCAM3 on the UK In- 
frared Telescope (UKIRT) in the K-band for approximately 
54 minutes, and in 20 of the 28 cases also in the J-band. 

For 5 sources (3C13, 3C41, 3C49, 3C65, 3C340) a fur- 
ther 3 to 4 hours of K-band observations have subsequently 
been taken using IRCAM3 in September 1998 (Best et al., 
in preparation), and these were combined with the origi- 
nal data to provide much deeper images. These further data 
were taken in the same manner as the original IRCAM3 data 
except for the use of the tip-tilt system that was available 
on UKIRT then, but wasn't available for the original runs; 
use of this system provided a significant reduction in the 
effective seeing. 

The HST data were reduced according to the stan- 
dard Spa ce Telescop e Science Institute (STScI) calibration 
pipeli ne (Lauer 1989| ), following the description of Best et al. 
( [1997| ). The UKIRT data were reduced using IRAF and 
mosaicked using the dimsum^ package, again following the 
general procedure outlined by Best et al. (1997). During 
the mosaicking process, the individual galaxy frames were 
block replicated by a factor of 4 in each dimension to allow 
more accurate alignment. The final images were then block- 
averaged in 2 by 2 pixels, resulting in frames with a pixel 
size of 0.142 arcseconds. The seeing varied from image to im- 
age, from about 0.7 to 1.2 arcsec. dimsum also produces an 
output exposure map indicating the exposure for each pixel 
on the combined image; this is useful for weighting the dif- 
ferent regions of the images, which reach different limiting 
sensitivities owing to the dithering process used in taking 
the data (see Best et al. 1997). 

This dataset has a number of advantages and disadvan- 
tages for use in studies of clustering around distant radio 
galaxies. On the negative side, since these data were origi- 
nally taken with the goal of studying the host galaxies of the 
radio sources, no comparison fields were taken with which 
the data can be compared to statistically remove background 
counts. Further, IRCAM3 is only a 256 by 256 array with 
a field of view per frame of just over 70 by 70 arcseconds. 
After mosaicking this provided a final image of nearly 100 
by 100 arcsec, which corresponds to a little over 800 kpc at 
a redshift of one for the adopted cosmology, but the high- 
est sensitivity is only obtained in the central 50 by 50 arc- 
sec region. This is a significantly smaller size than a cluster 
would be expected to have at this redshift, and so only those 
galaxies in the central regions of any prospective clusters will 
be investigated (cf. the much larger extent of the structure 
around the radio galaxy 3C324 at redshift 1.206, as demon- 
strated in Figure 2 of Dickinson (1997), although just over 
half of the associated cluster members found there would lie 
within the IRCAM3 field of view). 

Although these two factors limit the usefulness of the 
dataset, they are far outweighed by the benefits. First, the 



i DIMSUM is the 'Deep Infrared Mosaicking Software' package de- 
veloped by Eisenhardt, Dickinson, Stanford and Ward. 



large size and essentially complete nature of the sample 
make it ideal for surveying the average environments of 
these sources. Second, the multi-colour nature of the dataset 
provides important information on any cluster membership 
through the creation of colour-magnitude diagrams and the 
determination of photometric redshifts. Third, the availabil- 
ity of near-infrared data is essential for such studies at these 
high redshifts since it continues to sample the old stellar pop- 
ulations longwards of the 4000A break; in the K-waveband 
the K and evolutionary corrections are small and relatively 
independent of morphological type, even at redshifts z ^ 1, 
in contrast to the situation at optical wavelengths. Fourth, 
the availability of the high resolution HST data allows accu- 
rate star-galaxy separation down to the faintest magnitudes 
studied. For these reasons, this dataset is well-suited to in- 
vestigating the environments of z ~ 1 radio galaxies. 

2.2 K— band image detection and photometry 

Throughout these analyses, the K-band data were used as 
the primary dataset. Image detection and photometry on the 
K-b and frames was carried o ut using SEXTRACTOR version 
2.1 (Berlin & Arnouts 1996). As a result of the dithering 
technique used to obtain near-infrared data, the exposure 
time varies with position across the image meaning that the 
rms background noise level also varies. SEXTRACTOR allows 
the supply of an input weight map by which the local detec- 
tion threshold is adjusted as a function of position across the 
image to compensate the varying noise levels, thus avoiding 
missing objects in the most sensitive regions of the image or 
detecting large numbers of spurious features in the noisiest 
regions. The output exposure map of each field produced 
by DIMSUM was used as the weight map for SEXTRACTOR. 
The source extraction parameters were set such that, to be 
detected, an object must have a flux in excess of 1.5 times 
the local background noise level over at least A'^ connected 
pixels, where A'^ was varied a little according to the seeing 
conditions, but in general was about 20 (equivalent to 5 pix- 
els prior to the block replication and averaging during the 
mosaicking procedure). 

To test the validity of this extraction method, a search 
for negative holes was carried out using the same extrac- 
tion parameters, and resulted in a total of only 43 negative 
detections throughout the 28 images, that is, about 1.5 neg- 
ative detections per frame. The fluxes associated with these 
negative detections all correspond to positive features below 
the 50% completeness limit and so it is expected that there 
are essentially no false positive detections above that limit. 

The output from SEXTRACTOR was examined carefully, 
with minor (typically 1-2 per field) modifications being 
made to the catalogue if necessary: objects coinciding with 
the spikes of bright stars were removed, the data for occa- 
sional entries that had been separated into two objects by 
SEXTRACTOR but were clearly (on comparison with data at 
other wavebands) a single object were combined, and very 
occasionally it was necessary to add to the catalogue an ob- 
vious object which had been missed due to its proximity 
to a brighter object. Further, all objects that lay within 21 
pixels (3 arcsec) of the edge of an image were removed from 
the catalogue in case their magnitudes were corrupted by 
aperture truncation. SEXTRACTOR provides a FLAG param- 
eter to indicate the reliability of its measured magnitudes. 
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After this trimming of objects close to the edges of the fields, 
less than 2% of all objects had high flag values indicative 
of truncated apertures, saturated pixels, or corrupted data. 
These objects were also removed from the catalogue. 

SEXTRACTOR's mag_best estimator was used to deter- 
mine the magnitudes of the sources; thi s yiel ds an estimate 
for the 'total' magnitude using Kron's (1980) first-moment 



algorithm, except if there is a nearby companion which may 
bias the total magnitude estimate by more than 10% in 
which case a corrected isophotal magnitude is used instead. 
The determined magnitudes were also corrected for galactic 
extinc tion, using the extinction maps of Burstein and Heiles 
( [l982| ). 

To investigate the accuracy of these total magnitudes, 
and the completeness level of the source extraction as a func- 
tion of position on the image, Monte-Carlo simulations were 
carried out using the following 4-step process. 
Step 1: The point-spread function (PSF) of each K-band 
image was determined using objects which were unsaturated 
and which were unresolved on the HST images. 
Step 2: A series of model galaxies was made. 30% of 
the galaxies were assumed to be ellipticals, with radial 
light profiles governed by de Vaucouleurs' law, I[r) oc 

2 , /, /, N2\l/8l 



exp 



[{a/a^Y + {h/b^yy'" , where a and b are the 

distances along the projected major and minor axes, and Oc 
and 6c are the characteristic scal e leng ths in those directions. 
Lambas, Maddox and Loveday (1992) provide a parameter- 
ization for the ellipticities of elliptical galaxies in the APM 
Bright Galaxy Survey in terms of the parameter p — fec/tSc. 
An ellipticity was drawn from each galaxy at random from 
this distribution, and the position angle of the major axis 
was also chosen randomly. Songaila et al. ( 1994 ) found that 



for galaxies with 18 < /i' < 19 the median redshift is about 
0.6, which means that for K > 18 (where detection and 
completeness are being tested) the majority of the galax- 
ies will be distant enough that their angular sizes are rel- 
atively insensitive to redshift. High redshift ellipticals have 
typical characteristic sizes ranging from 2 to lOkpc (e.g. 
Dickinson 1997), and so the apparent characteristic size, rc 
{— (ocfoc)^^^), of each elliptical was chosen randomly from 
the range 0.2 to 1.2 arcsec. 

The remaining 70% of the galaxies were built using 
a single exponential profile appropriate for galaxy disks, 
7(r) oc exp[— r/rd], where rd is t he cha racteristic scale length 
of the disk. Mao, Mo and White (1998) show that the charac- 
teristic scale lengths of disks decreases w ith re dshift approx- 
imately as [1 + z)~^ , but Simard et al. (199!:) shows that a 
magnitude-size correlation somewhat counter-balances this: 
at a given apparent magnitude, higher redshift objects must 
be intrinsically brighter and so are larger. From their data, 
the typical scale length of galaxy disks at around our com- 
pleteness limit will be 1-4 kpc. The disk scale length of each 
galaxy was therefore chosen at random from the range 0.1 to 
0.5 arcsec. The disk inclination was chosen at random from 
to 90 degrees, with inclinations greater than 75 degrees being 
replaced by 75 degrees to account for the non-zero thickness 
of the disk. The projected orientation of the galaxy was also 
chosen randomly. 

The model galaxies were then convolved with the stellar 
PSF derived in step 1. 

Step 3: Five stellar or galactic objects were added to each 




Figure 1. The completeness fraction for both stellar (dashed 
line) and galactic (solid line) model objects as a function of input 
magnitude, averaged over all 28 fields. 



frame with a random position and scaled to a random mag- 
nitude in the range 15 < K < 22. SEXTRACTOR was then 
run on the new image using the same input parameters as 
for the original source extraction, to see if the added objects 
were detected, and if so to determine the difference between 
the total magnitude measured by SEXTRACTOR and the true 
total magnitude. 

Step 4- Step 3 was repeated until 25000 stellar objects and 
25000 model galaxies had been added to each image. 

From these results, the mean completeness fraction was 
determined as a function of magnitude for both stars and 
galaxies over the entire set of images. The results are shown 
in Figure 0. The 50% completeness limits are K = 19.75 for 



galaxies and K = 20.4 for stars. As discussed in Section 2.1 



new deeper observations are available for five fields, and con- 
sidering these deeper fields alone, the 50% completeness lim- 
its are 20.2 and 20.9 for galaxies and stars respectively. 

The input and measured magnitudes of each detected 
object can also be compared. Typical results for an individ- 
ual field are shown in Figure |2| Note the generally low scat- 
ter, particularly at bright magnitudes, but the occasional 
source with a large deviation. These are caused by the prox- 
imity of the model object to a bright source in the field; 
large errors such as this are avoided in the observed data by 
careful examination of the output SEXTRACTOR catalogues, 
as discussed above. The mean difference between the input 
and measured magnitudes, and the scatter of the measured 
magnitudes around this mean, are shown as functions of in- 
put magnitude for both model stars and galaxies in Figure ^ 

2.3 Multi-colour source extraction 

The J-band frames were aligned with the K-band frames by 
using a number of objects that appeared unresolved on the 
IfST images. This alignment often involved a small shift in 
position of the frames, but generally much less than a degree 
of frame rotation. SEXTRACTOR was then run in its double 
image mode using the J and K-band images. In double im- 
age mode SEXTRACTOR uses one image to detect objects and 
define the apertures to be used for flux determination, and 
then the fluxes and magnitudes are measured from the sec- 
ond image. In this way J-band magnitudes or upper limits 
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Figure 2. The distribution of measured magnitudes as a function 
of input magnitudes for the model stars (upper plot) and galaxies 
(lower plot). The plots are for 1000 model objects. No point is 
plotted if the model object is not detected by SEXTRACTOR. 



Figure 3. The mean (upper plot) and rms scatter (lower plot) 
of the difference between the measured magnitudes and input 
magnitudes for the model stars (dashed line) and galaxies (solid 
line). 



were determined for the sources detected in the K-band, 
through exactly the same apertures. 

For the HST data, SEXTRACTOR was first run on the 
HST frames obtained at the end of the calibration proce- 
dure. Then, the four separate WFPC2 frames were overlaid 
individually with the K-band data, using between 3 and 15 
unresolved objects visible at both wavelengths. In this pro- 
cess the HST data was re-pixelated to match the UKIRT 
data, essential for a combined running of SEXTRACTOR on 
the K-band and HST images; the re-pixelated and aligned 
HST images were then convolved to the angular resolution of 
the K-band data using a Gaussian convolving function. SEX- 
TRACTOR was run in double image mode on the convolved 
HST data, thus measuring accurate fluxes and magnitudes 
for the objects through the same apertures as the K-band 
data. The HST fluxes were corrected for the small differ- 
ences in gain ratios between the different WFPC2 chips and 
using a 4% linea r correction ramp fo r the charge transfer 
efficiency effects (Holtzman et al. 1995). They were also cor- 



rected for galactic e xtinct ion using the extinction maps of 
Burstein and Heiles (1982). Finally, for each object the 'stel- 



laricity index' (see Section 2.4), used for the separation of 



stars and galaxies, was replaced by that calculated in the run 
of SEXTRACTOR on the original HST frame: using the high- 
est angular resolution data provides a much more accurate 
determination of this parameter. 



For both the J-band and HST frames, only objects 
which were detected with fluxes of at least 3a were con- 
sidered, where a is the uncertainty on the flux measurement 
provided by SEXTRACTOR; this flux error estimate includes 
the uncertainty due to the Poisson nature of the detected 
counts and that from the standard deviation of the back- 
ground counts. An additional source of flux error arises from 
the uncertainty in the subtraction of the background count 
level as a function of position across the image. This value 
was estimated as the product of the area of the extraction 
aperture and the rms variation of the subtracted background 
flux across the image. This background subtraction error es- 
timate was combined in quadrature with the flux error given 
by SEXTRACTOR to determine the uncertainties on the mag- 
nitudes of the extracted objects. 



2.4 Star— galaxy separation 

SEXTRACTOR. provides a 'stellaricity index' for each object, 
which is an indication of the likelihood of an object being 
a galaxy or a star, base d on a neural network technique 
(Bertin & Arnouts 1996). In the ideal cases a galaxy has 
a stellaricity index of 0.0 and a star has 1.0. In practice. 



at faint magnitudes, low signal-to-noise and galaxy sizes 
smaller than the seeing lead to an overlap in the calculated 
stellaricity indices for the two types of object. 
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Figure 4. A plot of stoUaricity index against J — K colour: stars 
generally have bluer colours and higher stellaricity indices than 
galaxies, and so combining these two parameters allows fairly ac- 
curate star— galaxy separation. The stars lie to the upper left of 
the line. The field plotted here is 3C22, chosen for this demon- 
stration because of its high fraction of stars. 



For sources which were present on the HST frames (typ- 
ically only about 80% of objects, the precise percentage de- 
pending upon the sky rotation angle at which the HST frame 
was taken) , the stellaricity index for each object determined 
from the un-convolved HST data was adopted: the high res- 
olution of these images meant that, apart from the very red- 
dest objects, star-galaxy separation|^ was relatively unam- 
biguous right down to the completeness limit of the K-band 
frames. For the objects for which no HST data was avail- 
able, stellaricity indices were taken from the K-band data. 
At magnitudes K < 17.5, star-galaxy separation could be 
carried out directly from these stellaricity values. At fainter 
magnitudes it was still possible to provide a fairly accurate 
segregation of galaxies from stars by combining the stellar- 
icity index with the J — K (or F814W— /C) colour of the 
object, as shown in Figure W: the bluer average colour of 
stars can be used to separate out the stars and galaxies in 
the ambiguous range of stellaricity indices from about 0.5 to 
0.9. 

Figure ^ shows the total star counts as a function of 
magnitude derived from all of the frames, fit with the func- 
tion \ogN{K) = 0.247iC - 0.830. This good fit to the 
differential star counts using a single power-law distribu- 
tion clearly demonstrates that the star-galaxy separation is 
working well. 



3 GALAXY COUNTS 

Using the differences between true and measured magni 



2.2 



tudes determined for the simulated galaxies in Section 
(and Figure 3), the measured 'total' magnitudes of the 
galaxies were converted to true 'total' magnitudes, and these 
were binned in 0.5 magnitude bins over the magnitude range 
14 < < 20 to determine the raw galaxy counts, riraw. 



S Of course, star— galaxy separation using the resolved or unre- 
solved nature of sources means that any quasars will be placed 
into the star category, but their contribution is negligible. 



Figure 5. The differential star number counts in the K-band, 
in AK = 0.5 magnitude intervals, fit using a single power law 
distribution. 



These were corrected for completeness using the observed 
mean completeness fraction as a function of magnitude from 
the simulations (Figure |l]) . A possible source of systematic 
error in the number counts at faint magnitudes must also be 
taken into account: the increased scatter in the photometric 
measurements at faint magnitudes, and the fact that there 
are more faint galaxies than bright galaxies, mean that it 
is more likely for faint galaxies to be brightened above the 
K — 20 limit (or be brightened and move up one magnitude 
bin) than for brighter galaxies to appear erroneously in a 
fainter bin or fall out of the catalogue completely. This leads 
to an apparent increase in the number counts at the faintest 
magnitudes. The results from the model galaxy simulations 
were used to correct for this effect: the required corrections 
were smaller than 10% in all cases. Combining these effects 
produced corrected galaxy counts (ric), and scaling by the 
observed sky area, produced final counts per magnitude per 
square degree, A''c. 

Number counts were also derived using the same tech- 
nique c onsi dering just the five fields with deeper images (see 
Section 2.1). This provided determinations of galaxy counts 
for 17.5 < K < 20.5 which may be more reliable at the 
fainter magnitudes due to smaller completeness corrections. 

These galaxy counts are tabulated in Table |^ and plot- 
ted in Figure ^. The error on the number counts in each bin 
was calculated from two factors: the Poissonian error on the 
raw galaxy counts, and an error in the completeness correc- 
tion which is generously assumed to be 30% of the number 
of counts added in the correction. The number counts are 
compared on the plot with counts from various K-band field 
surveys in the literature. The final column in Table ^ shows 
the mean galaxy counts per magnitude bin determined from 
these literature field surveys. 

The galaxy counts derived in this paper show an excess 
of counts relative to the literature counts, which for magni- 
tudes K > 15.5 is at greater than the la significance level. 
For comparison, the K-magnitudes of the radio galaxies 
themselves span the range 15.4 < K < 18.0. The integrated 
excess counts over the magnitude range 15.5 < K < 20 cor- 
responds to, on average, 11 galaxies per field. Unfortunately, 
as described in Section Ejl, the original goals of these obser- 
vations was to carry out detailed studies of the radio source 
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Table 1. K— band galaxy counts as a function of magnitude. 
The columns give the raw counts (nraw ) , the counts corrected for 
completeness and magnitude biasing effects (ric), the corrected 
counts per square degree per unit magnitude (Nc), the error on 
this, and the mean literature counts per souare degree per unit 
magnitude from the data shown in Figure M (Afm). 
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Figure 6. The K— band galaxy number counts per magnitude 
per square degree, in AK = 0.5 magnitude intervals. The solid 
symbols represent the data presented in this paper and tabu- 
lated in Table ^ for X < 19 the values plotted are those from 
'all fields' and for K > 19 they are those from the 5 deep fields 
only. The other symbols represent data from the literature: open 
squares and (point-up) triangles are respectively the bright and 
faint samples of Minezaki et al. (1998); diamonds are from Djor- 
govski et al. (1995); circles are from Moustakas et al. (1997); 
x's are from Szokoly et al. (1998); asterisks, inverted triangles 
and stars are respectively the Hawaii Medium Deep Survey, the 
Hawaii Deep Survey and the Hawaii Medium Wide Survey pre- 
sented by Gardner, Cowie and Wainscoat (1993); open crosses 
are from McLeod et al. (1995); bow-ties (m) are from Bershady, 
Lowenthal and Koo (1998). The solid line repr esen ts a passive 
luminosity evolution model, discussed in Section 4.3. 



host galaxies, and so no blank sky frames were taken. There- 
fore no comparison can be made between radio galaxy and 
blank sky frames to determine whether the excess counts 
found here are associated with structures surrounding the 
radio galaxies or whether there is some systematic offset be- 
tween the counts calculated here and the literature counts. 

We believe that the majority of the excess counts are as- 
sociated with the presence of the radio galaxies, for a n umber 
of reasons. First, Roche, Eales and Hippelein (1998) found 
a similar excess in the K-band galaxy counts in a study 
of the fields surrounding 6C radio galaxies at similar red- 
shifts. Second, as reviewed in the introduction, there are 
numerous lines of evidence suggesting that there are clus- 
ters around at least some of these objects. Third, later in 
this paper it will be shown that cross-correlation analyses 
and colour-magnitude relations indicate an overdensity of 
galaxies comparable to that observed. 



4 CROSS-CORRELATION ANALYSES 

4.1 Angular cross-correlation estimators 

The clustering of galaxies around the radio galaxies can 
be investigated using the angular cross-correlation function 
1^(9), which is defined from the probability of finding two 
sources in areas SQi and 5i}2 separated by a distance 6: 

sp = N^[i + w{e)]5niSn2 

where TV is the mean surface density of sources on the sky. 
From this, it follows that for a given survey: 

DZ?(e) = ino(no-l)[l+™(e)]^^^^, (1) 

where DD{9) is the number of data-data pairs with angular 
separation of between 9 and 6 + 59, no is the total number 
of sources in data catalogue, SI is the total angular area 
of sky sampled and {Sfl-o{9)) is the mean angular area of 
sky accessible at a distance 9 to 9 + S9 around the data 
points. Clearly {Si}o{9)) is extremely difficult to calculate 
due to boundary effects; various estimators for w{9) have 
therefore been derived using comparisons between the data 
points and catalogues of randomly distributed points (e.g. 
see Cress et al. 1996 for a discussion). For the estimator 
adopted in this paper, it is considered that if na random 
points are added to the image (where na ^ no to minimise 
errors introduced by the random catalogue) then the number 
DR{9) of data-random pairs between 9 and 9 + 59 will be 
given by: 

{5no{9)) 



DR{9) = nonn- 



(2) 



Combining equations ^ and ^ gives the following estimator 
for w{9): 



2nn DD{9) 
{no - 1) DR(9) 



(3) 



In addition, it is possible to consider the angular cross- 
correlation function just around an individual source, in 
this case the radio galaxy. For this, the parameters be- 
come RGD{9) = {no ~ + Wr!,{9)]{5nnG{9)} /il and 
RGR{9) = nR{5QRG{9))/n, where RGD{9) and RGR{9) 
are respectively the number of radio galaxy-data and radio 
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galaxy-random pairs between and 8 + 59, Wrg{0) is the 
angular cross-correlation function for galaxies around the 
radio galaxy, and {SnRG{0)} is the angular area of sky ac- 
cessible at a distance 6 to 9 + 66 from the radio galaxy. These 
equations combine to give: 



_ nR RGDje) 
' no RGR(e) 



(4) 



4.2 Calculating w(6l) and Wrg{0) 

One problem with estimating ■w{9) and Wrg(6) from the cur- 
rent data is that, due to the K-band dithering technique, 
the detection limit is a function of position across each im- 
age. More objects are detected in the central regions of the 
image, which would produce a spurious peak in the angu- 
lar cross-correlation statistics. This non-uniform noise level 
can be accounted for using the simulations which were dis- 
cussed in Section 2.2. From these simulations, for each field 



the completeness fraction can be calculated as a function 
of both magnitude and position, and then applied to the 
random catalogue in the following way. 
Step 1: Excluding those objects determined to be stars, 
DD{6), RGD{6), and n-o were determined for each frame 
down to a chosen limiting magnitude ifiim, for a set of bins 
in e. 

Step 2: 25000 random objects were each assigned a random 
position on the image, and a magnitude K < Kun^ drawn at 
random from a distribution matching the observed number 
counts (Figure ^). 

Step 3: For each random object, the completeness fraction 
of objects of that magnitude and position (averaged in 50 
by 50 pixel bins) was determined from the simulations. The 
object was accepted or rejected at random with a likelihood 
of acceptance based upon the derived completeness fraction. 
Step 4- For those objects which remained in the catalogue, 
DR{6), RGR{6), and jir were calculated. Hence w{6) and 
WTg{9) were derived. 

This process was repeated for all images and for two 
different limiting magnitudes, ifiim = 19 and Kum ~ 20. 
Kum = 19 corresponds to the limit at which the data are still 
almost 100% complete, whilst Ku^n = 20 requires a signifi- 
cant completeness correction but provides improved statis- 
tics for the galaxy counts. The results were averaged over 
all 28 fields; investigating the values for individual fields was 
not attempted since for these small fields of view the mea- 
sured amplitudes of single images are too strongly affected 
by variations in the background (and foreground) counts to 
provide results with a high statistical significance; only by 
combining the data of the 28 fields is a sufficiently robust 
measurement obtained. An estimation for the uncertainty 
in the value of the combined w{6) in each bin was provided 
by the scatter in the values between the different fields. The 
results are shown in Figure ^. 

w{6) is usually assumed to have a power-law form; 
w(6) — A{9/deg)~^ . If so, the observed w(6) will follow a 
form w{6) — A{6/deg)~^ — G, where G is known as the inte- 
gral constraint and arises from the finite size of the field 
of view. Its value can be estimated by integrating w{9) 
over the area of each field, and for our data corresponds 
to Cgg = 41.1j4gg and Crg = 48.0j4rg for, respectively, all 
galaxy-galaxy pairs and just radio galaxy-galaxy pairs. It is 



Table 2. Amplitudes of the angular and spatial cross— correlation 
functions for the fields surrounding the SCR galaxies. Both the 
amplitudes for all galaxy-galaxy pairs and those considering only 
radio galaxy-galaxy pairs are given, down to two limiting magni- 
tudes. For radio galaxy-galaxy pairs these are convert ed t o spatial 
cross— correlation amplitudes as described in Section 
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Sources 




A 


AA 


B 


AB 


All g-g pairs 


K <19 


0.0031 


0.0012 








K <20 


0.0015 


0.0004 






Only rg-g pairs 


K <19 


0.0092 


0.0037 


600 


240 




K <20 


0.0093 


0.0021 


510 


120 



not possible to determine both the amplitude and the slope 
of the fit from the current data, and so the canonical value of 
5 = 0.8 (which also seems to be appropriate at high redshifts; 
Giavalisco et al. 1998) has been adopted. The observed data 
in Figure ^ has been fit with functions of this form, and the 
resulting angular cross-correlation amplitudes are provided 
in Table H 

The results at if < 19 and if < 20 are in approxi- 
mate agreement. For the galaxy-galaxy pairs there may be 
a decrease of Agg as fainter magnitude limits are used (as 
expected, e.g. see Roche et al. 1998), but the errors are too 
large to determine this with any degree of confidence. The 
galaxy-galaxy amplitudes are similar to those derived by 
Roche et al. to the same magnitude limits in the fields of 
radio galaxies at z ~ 0.75. A more obvious feature is that 
at both magnitude limits the cross-correlation amplitude 
around the radio galaxy is significantly larger (Arg S> Agg). 
Figure ^ shows that much signal originates at small angu- 
lar separations {9 < 10") and so this cannot be related to 
any problems with completeness in the outer regions of the 
frames; the similarity of the results for is' < 19 and A" < 20 
also demonstrates this. 



4.3 The spatial cross— correlation amplitude 

As has been described by many authors (e.g. Longair & 
Seldner 1979, Prestage and Peacock 1988) it is possible to 
convert from an angular cross-correlation amplitude to a 
spatial cross-correlation amplitude if the galaxy luminosity 
function is known. The spatial cross-correlation function is 
usually assumed to have a power-law form: 



( — 

V Mpc 



where the power-law slope is related to the slope of the 
angular cross-correlation function by 7 = 5 -I- 1 (=1.8 for 
the 5 = 0.8 adopted here). The spatial cross-correlation 
amplitude, i3rg, is then related to that of the angular cross- 
correlation by Aj-g = H{z)B^g (see Longair & Seldner (197t ) 
for a full derivation) where: 



H{z) 



Ng 



Djz) 
1 + z 



3-7 



0(mo,z) (7r/180) 



-(7-1) 



(5) 



Here, I-y is a definite integral which for 7 = 1.8 has a value 
of 3.8, A'g is the measured sky density of objects above the 
magnitude limit (per steradian), D{z) is the proper dis- 
tance to a source at redshift z, <f){mo,z) is the number of 
galaxies per unit comoving volume which at redshift z are 
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Figure 7. The angular cross-correlation function for all galaxy-galaxy pairs (left) and for just the radio galaxy-galaxy pairs (right) 
averaged over all of the UKIRT K— band images. The upper plots include all galaxies with K < 19, to which level the observations are 
almost complete; the lower plots include galaxies down to A" < 20, at which level there is a significant completeness correction, but the 
number statistics are increased. The plots are fit with functions of the form w{9) = A9~'^'^ — C, where C is the integral constraint and 
has a value of Cgg = 41.1Agg for all galaxy— galaxy pairs, and Crg = 48.0Arg for just radio galaxy-galaxy pairs (see text). 



more luminous than apparent magnitude mo, and the fac- 
tor (7r/180)~'^~^' is required to convert Arg from degrees to 
radians. The value of H{z) is only a very weak funct ion of 
7 provided 7 is of order 2 (Prestage & Peacock 1988), and 
so the fixing of 7 at 1.8 will introduce no significant errors. 

0(mo, 2) requires assumptions for both the K-correction 
of galaxies and the local luminosity function, as well as being 
cosmology dependent, and so this conversion wil l alwa ys be 
somewhat uncertain. In this paper a Schechter (1976) form 



is adopted for the local luminosity function, that is: 

^(L)d.^,-(A)%.p(_^)a(ii) 

where (j>{L)dL is the number of galaxies per comoving cubic 
Mpc with luminosity L to L + dL, (p* is the density normal- 
ization factor, L* is the characteristic luminosity and a is the 
faint-end slope. A pure luminosity evolution model is also 
adopted; that is, values for (j>* , 01 and L*(z = 0) are taken, 
and the only evolution of this function is then evolution 
of the value of L* with redshift in accordance with passive 
evolution predictions. As discussed in the introduction, pas- 
sive evolution models provide a good description of cluster 



galaxy properties back to 2 ~ 1, and they can also provide a 
reasonable fit to the observed K-band field number counts 
(see below). These models are, however, undoubtedly a sim- 
plification in view of hierarchical galaxy formation theories. 
In hierarchical models, fewer luminous galaxies are predicted 
to exist at high redshifts (e.g. see Figure 4 of Kauffmann and 
Chariot 1998), resulting in the values of (j){mo, z) and conse- 
quently H{z) being lower, and hence Brg will be increased. 
The pure passive luminosity approach therefore provides a 
conservative lower estimate for the spatial cross-correlation 
amplitude around high redshift radio galaxies. 

For the passive evolutionary models, the galaxy popula- 
tion was split into four different galaxy types: ellipticals and 
SO's (E's), Sa and Sb types (Sab's), Sc types (Sc's) and Sd 
types and irregulars (Sdm's). Galaxies of these types were 
built up using the Bruzual and Chariot (1993) stellar syn- 



thesis codes (1996 version), assuming a Scalo initial mass 
function, solar metallicity, a formation redshift zf = 10, and 
four different star formation histories (cf. Guiderdoni and 
Rocca-Volmerange 1990). The E's were assumed to form 
their stars in a rapid early burst with the star formation 
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decreasing exponentially on a rgfr = 0.5 Gyr timescale. The 
Sab's had 50% of their stars in a bulge component formed 
in the same manner as the E's, and the remaining 50% 
in a disk-like component with a much longer star forma- 
tion timescale {ts^ = 6 Gyr). The Sc's were modelled using 
just a single long star formation timescale {tsIt = 8 Gyr), 
and the Sdm's using a still longer timescale (rgfr — 50 Gyr) 
to produce significant star formation at the current epoch. 
These star formation histories approximately reproduce the 
colours of the various morphological types at the current 
epoch. It was assumed that all of these four morphologi- 
cal types displayed the same luminosity function (that is, 
M* and a), with the weightings of the different types (ie. 
the relative contributions to (p*) being 30%, 30%, 30% and 
10% respectively. In fact, at wavelengths as long as the K- 
band, the K and evolutionary corrections are small and rel- 
atively independent of morphological type, and so the re- 
sults are not strongly dependent upon any of these assump- 
tions. Changing the values adopted above by a sufficiently 
large amount that they produce unacceptable distributions 
of galaxy colours at the current epoch produces only 10-20% 
changes in the final K-band luminosity function. 

Local K-band lumi nosity functions have been derived 
by Ga rdner et al. (1993) , by Mobasher, Sharpies and Ellis 
(1993) and by Loveday (200C). All three of these datasets 
are approximately consistent with a characteristic absolute 
magnitude Mj^ = —25.1 and a slope a = —1.0. To inves- 
tigate whether the combination of this luminosity function 
and the assumption of passive evolution models produces 
acceptable results, and to calculate the appropriate value 
of (p* , the K-band number counts as a function of appar- 
ent magnitude that would be expected in this model were 
derived. The value of (j)* was adjusted to provide a good 
fit to the observed data points. The resulting fit, for a value 
(j)* ~ 0.004, is shown as the solid line in Figure^, and demon- 
strates that this simple model can provide a reasonable fit 
to the observed number counts. 

With this model luminosity function, it is possible to 
calculate (j}{mo,z) and, using the values of A'^g from Ta- 
ble (iVg = 9.2 X 10'^ sr"^ for K < 19 and 2.0 x 10*sr"^ 
for K < 20), hence to determine H{z). For K < 19, H{z) 
ranges from 15.4 x lO"'^ to 1.8 x 10"^ over the redshift range 
0.6 to 1.8. The mean value of 1/H{z) averaged over all 28 
radio source redshifts is 1/H{z) ^ 1.7 x 10®, and using this 
value to convert from A^g to Brg gives Big = 600 ± 240. For 
K < 20 the corresponding range of H{z) is from 12.2 x 10^^ 
to 3.1 X 10"^, with a mean value of 1/H{z) ^ 1.4 x 10®, 
corresponding to Brg = 510 ± 120. 

These values can be interpreted physically by com- 
paring with the equivalent values for Abell clusters cal- 
culated between the central galaxy and the surrounding 
galaxies (-Beg). This has been calculated independently as 
a function of Abell cluster r ichness by a number of authors 
( iPrestage fc Peacock 19881 |Hill fc Lilly 199l| ; [Andersen & 
Owen 1994; Yee fc Lopez-Cruz 1999| ). Converting their val 



ues to Ho — 50 km s ^ Mpc ^ and 7 = 1.8 they are all 
approximately consistent with each other, and average to 
Scg fs 350 for Abell class and Beg « 710 for Abell class 1. 
The environments surrounding the redshift one radio galax- 
ies are therefore comparable, on average, to those of clusters 
of between Abell Classes and 1 richness. 



Hill and Lilly (1991) further showed that there is a cor- 
relation between the value of Beg and the parameter A^o.s, 
where is an Abell-type measurement defined as the 

net excess number of galaxies within a radius of 0.5 Mpc 
of the central galaxy with magnitudes between mi and 
mi -I- 3, mi being the magnitude of the central galaxy. Us- 
ing a larger dataset, Wold et al. (2000) calibrated this re- 
lation as Beg = (37.8 ± 10.9)Ai'o.5. The average Brg value 
for K < 20 then implies that an average net excess of 13.5 
galaxies should be found around each radio galaxy within a 
radius of 0.5 Mpc and with a magnitude down to three mag- 
nitudes fainter than the radio galaxy. The data presented in 
this paper cover about 80% of this sky area and, since the ra- 
dio galaxies have typical magnitudes of -K" ~ 17, the galaxy 
counts to if = 20 do sample approximately 3 magnitudes 
below the radio galaxy. Therefore, the net excess counts of 
11 galaxies per field down to A' = 20 (see Section ^ is fully 
consistent with the determined value of Brg. 



5 COLOUR-MAGNITUDE AND 
COLOUR-COLOUR RELATIONS 

5.1 Near— infrared colour— magnitude relations 

For each field the magnitudes of all galaxies were deter- 
mined through all of the different filters, as described in 
Section After correction for galactic reddening, these 
were used to construct colour-magnitude relations for each 
field. The near-infrared J — K vs K relation for each field (if 
J-band data were not taken, the longest wavelength HST 
filter available was used instead) are shown in Figure ^ in 
order of ascending redshift. The stars (open diamonds) and 
galaxies (asterisks) have been separated on these diagrams 
using the technique described in Section and the radio 
galaxy is plotted using a large symbol. The uncertainties 
on the measured colours are not shown to avoid cluttering 
up the diagrams; at bright magnitudes {K < 17) they are 
dominated by calibration uncertainties (< 0.1 mags), but 
increase to about 0.3 magnitudes by if ~ 20. 

The colour-magnitude relation for the Coma cluster has 
been redshifted and evolved according to the passive evolu- 
tionary models for elliptical galaxies described in the pre- 
vious section; this redshifted relation is over-plotted on the 
figures. Note that the radio galaxies generally have near- 
infrared colours similar to this theoretical line, indicating 
that these are old elliptical galaxies (e.g. Best et al. 1998b), 
but not in all cases: some of the radio galaxies appear redder 
because a heavily reddened nuclear component contributes 
to the K-band emission (e.g. 3C22, 3C41; Economou et al. 
1995, Rawlings et al. 1995, Best et al. 1998), whilst for oth- 
ers the bluer filter is significantly contaminated by excess 
optical-UV emission induced by the radio source (the 'align- 
ment effect', e.g. McCarthy et al. 1987). 

A number of features are apparent from these near- 
infrared colour-magnitude diagrams. Considering initially 
those radio galaxies with redshifts 2 < 0.9, many of the 
sample show reasonably convincing evidence for associated 
clusters; here it should be born in mind that for an Abell 
Class 1 cluster only 10 to 15 associated cluster galaxies 
are expected in the observed region of sky down to three 
magnitudes below the magnitude of the radio galaxy (see 
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Figure 8. Near— infrared colour-magnitude relations for the SCR radio galaxies. The fields and redshifts are labelled above each diagrams. 
Galaxies are represented by asterisks and stars by open diamonds. Upper limits are plotted as an arrow for a galaxy and as an arrow 
plus a diamond symbol for a star. The large plotted symbol represents the radio galaxy. The solid lines display the colour-magnitude 
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Figure 8. (continued) relations that would be obtained by simply passively evolving the colour-magnitude relation observed in the 
Coma cluster back in redshift, assuming an elliptical galaxy which formed its stars in a short time-span at z = 10. 
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Section [4.3| ). Some fields (e.g. 3C34, 3C337) clearly show 
at least this number (the 'background' counts at these K- 
magnitudes and colours are small, as can be seen by a com- 
parison with the colour-magnitude relations of the higher 
redshift sources), whilst other fields (e.g. 3C217) show few if 
any associated galaxies. There are clearly large source-to- 
source variations in environmental richness. 

Where a clear colour-magnitude relation is observed, 
the mean colour of this relation lies close to that calculated 
theoretically by just passively evolving the Coma colour- 
magnitude relations back in redshift, and the scatter around 
the colour-magnitude sequence is small. These results have 
been shown before for optical and X-ray selected clusters 
at redshifts out to 2; ~ 0.8 (e.g. Stanford et al. 1998). 
There are some small deviations in the colour of the ob- 
served sequence from the passive evolution relation (e.g. see 
3C265), but none much larger than a couple of tenths of 
a magnitude. The near-infrared colours of the galaxies in 
the radio galaxy fields are therefore consistent with those 
observed in other clusters at the same redshift, implying 
that the excess galaxy counts are associated with a struc- 
ture at the radio galaxy redshift. This is important because 
of previous suggestions that powerful distant radio galaxies 
may be systematically amplified by foreground lensing struc- 
tures (Le Fevre et al. 1987; Bem'tez, Martmez-Gonzales and 
Martin-Mirones 1997); were such structures to be present, 
they could account for both the excess K-band counts and 
the peak in the cross-correlation statistics, but the colour- 
magnitude relations argue against this. 

At higher redshifts the evidence for clear colour- 
magnitude relations is much poorer. This is mainly because 
the combination of the greatly increased contribution from 
field galaxies at these fainter magnitudes and the increased 
scatter in the colours due to photometric uncertainties, re- 
sults in any colour-magnitude sequence appearing much less 
prominent. The difficulty of selecting cluster candidates at 
these redshifts on the basis of a single colour can b e gauged 
by examining 3C324; for this field, Dickinson ( 1997 ) has con- 
firmed the presence of a poor cluster of galaxies, but this is 
barely apparent from its colour-magnitude relation. 

Although no prominent colour-magnitude sequences 
are seen in the higher redshift fields, there remains a net 
excess of K-band counts: if the fields of the radio galaxies 
are divided into two redshift bins, no significant differences 
in the faint galaxy counts are seen between the high and low 
redshift fields. A simple analysis also shows that the excess 
in the fields of the higher redshift radio galaxies is comprised 
of red {J — K> 1.75) galaxies. In Table M are given the mean 
number of galaxies per field with magnitudes 17 < < 20 
and colours J - K > 1.75 or 1.25 < J - K < 1.75, in the 
low and high redshift bins. There are more galaxies with 
bluer J — K colours in the fields of the lower redshift radio 
galaxies than in those at higher redshifts, due to the associ- 
ated cluster galaxies in the lower redshift fields which have 
colours 1.25 < J - K < 1.75 (cf. Figure ||). On the other 
hand, there are more galaxies with J — K > 1.75 colours in 
the high redshift than low redshift fields. The excess K-band 
galaxy counts in the high redshift field appear to be predom- 
inantly associated with red galaxies, with colours similar to 
those expected for old cluster galaxies at these redshifts. 
This again indicates that the excess number counts are as- 



Table 3 . The mean number of galaxies per field with magnitudes 
17 < K < 20 and blue/red colours, as a function of redshift. 



2 < 0.9 
2 > 0.9 



No. fields 



9 

11 



Mean number of galaxies 
1.25 < J - < 1.75 J-K>1.75 
19.6 ±2.5 8.1 ±1.7 

14.1 ±1.8 13.1 ±2.6 



sociated with a structure at the redshift of the radio galaxy 
rather than a foreground structure. 



5.2 Multi-colour relations 

In Figure ^ are plotted a complete set of colour-magnitude 
relations for six galaxies in the sample, chosen to be those 
which show amongst the best examples of near-infrared 
colour-magnitude relations for their redshiftj^ From these 
it is apparent that, although the scatter around the near- 
infrared colour-magnitude sequence remains small even at 
these high redshifts, colours that reach shortward of the 
rest-frame 4000A break show a dramatic increase in the 
scatter of the relation, (e.g. compare the various relations 
for 3C34); these colours can be strongly influenced by small 
amounts of recent or on-going star formation, indicating 
that this may be common in these high redshift clusters. 
Whether such star formation is in some way connected with 
the presence of a powerful radio source in these clusters can- 
not be distinguished from these data. 

To properly investigate the nature of the galaxies in 
these fields, all of the colour information must be used si- 
multaneously to derive photometric redshifts and investi- 
gate star formation activity. This is beyond the scope of 
the current paper but will be addressed later (Kodama & 
Best, in preparation). Here, the use of the all multi-colour 
information simulataneously is merely demonstrated in Fig- 
ure ^ through colour-colour plots for these six fields. For 
each field, the near-infrared J — K colour of each galaxy 
is plotted against its optical-infrared colour. These data are 
then compared against theoretical evolutionary tracks for 
the four different passively evolving galaxy models consid- 
ered in Section |j (E's, Sab's, Sc's, Sdm's). 

For the lower redshift radio galaxies there is clearly a 
large concentration of galaxies with colours very close to 
those of the model elliptical galaxy at the redshift of the ra- 
dio source. Further, there is a distribution of galaxies with 
colours between this and the colours of the model spiral 
galaxies at that redshift. In contrast, for the higher redshift 
radio galaxies, no strong concentration of galaxies is seen 
close to the elliptical galaxy prediction, and the number of 
cluster candidates lying between the locations of the ellipti- 
cal and spiral model galaxies is smaller than that found in 
the lower redshift cases. Clearly, despite the excess -counts 
and red J — K galaxies in the fields of these high redshift 
radio sources, to accurately investigate cluster membership 
requires the construction of photometric redshifts using sev- 
eral colours measured with high photometric accuracy. 



^ Equivalent plots for the other fields are available from the au- 
thor on request, but are not reproduced here to save space. 
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Figure 9. Colour— magnitude relations for the 3CR radio galaxies from the UKIRT and HST data. Symbols as in Figure |^. The redshifted 
Coma colour-magnitude relations are only plotted on figures for which both filters are longward of the rest-frame 4000A break, since 
shortward of that wavelength small quantities of recent or on— going star formation can have large eff'ects on the colours leading to only 
poorly defined sequences. 
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6 DISCUSSION 

In this paper a number of pieces of evidence for cluster- 
ing around distant SCR radio galaxies have been presented. 
These can be summarised as follows: 

• The K-band number counts show an overdensity of 
faint galaxies in the fields of the radio galaxies, with a mean 
value of 11 excess galaxies per field. 

• This excess is comparable to the galaxy overdensity ex- 
pected for a field of view of this size centred on a cluster of 
approximately Abell Class richness. 

• Cross-correlation analyses show a pronounced peak 
in the angular cross-correlation function around the radio 
galaxies. 

• Assuming that the galaxy luminosity function under- 
goes pure passive luminosity evolution with redshift, the cor- 
responding spatial cross-correlation amplitude lies between 
those determined for Abell Class and Abell Class 1 clus- 
ters. 

• The galaxies in the fields of most of the lower redshift 
radio galaxies in the sample show clear near-infrared colour- 
magnitude relations with only small scatter. The colours of 
these sequences are in agreement with those of other clus- 
ters at these redshifts, indicating that the excess number 
counts and cross-correlation peak are both associated with 
a structure at the redshift of the radio galaxy. 

• There is considerably more scatter in the relations in- 
volving shorter wavelength colours, suggesting low levels of 
recent or on-going star formation in many of the galaxies. 

• At higher redshifts the colour-magnitude relations are 
less prominent due to increased background contributions, 



but there is a clear excess of galaxies with very red infrared 
colours. 



These features all provide strong evidence that distant 
radio galaxies tend to reside in rich environments. The num- 
ber counts, the cross-correlation statistics, and the colour- 
magnitude relations all complement the previous results 
from X-ray imaging, narrow-band imaging, spectroscopic 
studies and radio polarisation studies discussed in Section |l| 
A coherent picture that most, but not all, high redshift ra- 
dio sources live in poor to medium richness clusters has now 
been built. 

Taking the results at face value, the environmental rich- 
ness around these z ~ 1 radio galaxies is higher than that 
around powe rful r adio galaxies &t z ^ .5 ca lculated by 



Yates et al. ( |l989| ) and by Hill and Lilly ( |l99l| ). This sug- 
gests that the increase between 2 = and z — 0.5 in the 
mean richness of the environments surrounding FRII radio 
galaxies, found by those authors, continues to higher red- 
shifts. 

Some notes of caution must be added to this conclu- 
sion. First, from the variations in richness of the colour- 
magnitude relation at any given redshift, (e.g. compare 
3C217 and 3C226), it is apparent that there is a wide spread 
in the density of the environments in these fields. Although 
most show some evidence of living in at least group environ- 
ments, not all powerful distant radio galaxies lie in clusters. 
Second, a simple visual comparison with the extremely rich 
high redshift cluster MS1054-03 (z = 0.83, van Dokkum 
1999) is sufficient to demonstrate that even at high red- 
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Figure 10. Colour-colour plots for fields of the six radio galaxies presented in Figure ^ The large triangle represents the radio galaxy, 
and the smaller upside-down triangles represent the other galaxies in the field. Galaxies with upper limits on their colo urs are indicated 
by arrows. The four tracks show the colour evolution of the passively evolving galaxy models described in Section \ E's - solid line; 
Sab's - dot-dash line; Sc's — dotted line; Sdm's - long dashed line. The tracks run from z = (lower left of plot) to 2: = 3, with the small 
crosses along each line representing steps of 0.1 in redshift, and redshifts 2 = 0, 1, 2, 3 being indicated by the larger crosses. 
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shifts, powerful radio galaxies still avoid the most extreme 
richness clusters. 



forming rich clusters. The signific ant scatter in the black 
hole versus bulge mass correlation (Kormendy & Richstone 



Fu ther, although the galaxy count excesses and the 1995) would, however, result in some scatter in the richness 



cross-correlation amplitudes have been compared with those 
of Abell clusters at low redshifts, in hierarchical galaxy for- 
mation models such comparisons will always be somewhat 
ambiguous. On-going mergers between galaxies mean that 
more sub-clumps are seen at higher redshifts, whilst the 
general galaxy cross-correlation length also evolves with red- 
shift in a manner dependent upon both cosmological param- 
eters and the method o f selecting the galaxy populations 
(Kauffmann et al. 199£). Therefore, quantitative interpre- 
tations of either of these parameters at high redshift must 
be considered with some care. On the other hand, in hier- 
archical growth models, the structures in which the radio 
galaxies reside will also continue to grow and evolve into 
much richer structures by a redshift of zero, meaning that 
the qualitative result that the high redshift radio galaxies 
lie in rich environments for their redshift is secure. 

Finally, it is important to consider the consequences for 
our understanding of the onset and nature of powerful ra- 
dio sources of a change in the environments of FRII radio 
galaxies from groups at low redsh ifts to clusters at high red- 
shift. As discussed by Best et al. (1998h), if this result holds 
then the standard interpretation of the tightness and slope 
of the Hubble K — z relation, that of 'closed-box passive 
evolution' of radio galaxies at 2: > 1 into radio galaxies at 
2; ~ 0, is no longer valid. It is not possible that the environ- 
ments can become less rich with progressing cosmic time. 
Instead, Best et al. propose that powerful radio galaxies se- 
lected at high and low redshift have different evolutionary 
histories but must contain a similar mass of stars, a few 
times 10^^ M0, so conspiring to produce the observed 'pas- 
sively evolving' K— z relation. In their model, powerful FRII 
radio sources are seldom formed in more massive galaxies 
(that is, in central cluster galaxies at low redshifts) because 
of the difficulties in supplying sufficient fueling gas to the 
black hole: in rich low redshift clusters the galaxies and gas 
have been virialised and take up equilibrium distributions 
within the cluster gravitational potential, the galaxies have 
high velocity dispersions greatly reducing the merger effi- 
ciency, and there is a dearth of gas-rich galaxies close to 
the centre of the clusters which might merge with, and fuel, 
the central galaxy. Thus, the formation of a powerful radio 
source in these environments is a rare event (but can still 
happen, e.g. Cygnus A). 

At high redshifts, radio galaxies can be found in (proto) 
cluster environments because these are not yet virialised, 
have frequent galaxy mergers, and have a plentiful supply 
of disturbed intracluster gas to fuel the central engine and 
confine the radio lobes. The central cluster galaxies will be 
amongst the most massive galaxies at these redshifts and so, 
from the correlation between black hole mass and bulge mass 
(e.g. Kormendy and Richstone 1995), will have the most 
massive black holes. The kinetic energy of the relativistic 
radio jets of distant SCR radio galaxies corresponds to the 
Eddington luminosity of a bla ck hole with M ~ 10^ — 10^ Mq 
(Rawlings & Saunders 1991), implying that these sources 



are fueled close to the Eddington limit. Therefore the most 
powerful radio sources will tend to be powered by the most 
massive central engines and hence be hosted by the most 
massive galaxies, which tend to be found at the centres of 



of the environments of the radio galaxies. The data presented 
in this paper are in full agreement with this model. 
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